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We consider the cosmology that results if our observable universe is a 3-brane in a higher di- 
mensional universe. In particular, we focus on the case where our 3-brane is located at the Z2 
symmetry fixed plane of a Z2 symmetric five-dimensional spacetime, as in the Hof ava-Witten model 
compactified on a Calabi-Yau manifold. As our first result, we find that there can be substantial 
modifications to the standard Friediuann-Robertson- Walker (FRW) cosmology; as a consequence, 
a large class of such models is observationally inconsistent. In particular, any relationship between 
the Hubble constant and the energy density on our brane is possible, including (but not only) 
FRW. Generically, due to the existence of the bulk and the boundary conditions on the orbifold 
fixed plane, the relationship is not FRW, and hence cosmological constraints coming from big bang 
nucleosynthesis, structure formation, and the age of the universe difficult to satisfy. We do wish 
to point out, however, that some specific choices for the bulk stress-energy tensor components do 
reproduce normal FRW cosmology on our brane, and we have constructed an explicit example. As 
our second result, for a broad class of models, we find a somewhat surprising fact: the stabilization 

jr \ ■ of the radius of the extra dimension and hence the four dimensional Planck mass requires unrealistic 

fine-tuning of the equation of state on our 3-brane. In the last third of the paper, we make remarks 
about causality and the horizon problem that apply to any theory in which the volume of the extra 

f^l . dimension determines the four-dimensional gravitational coupling. We point out that some of the 

^ ' assumptions that lead to the usual inflationary requirements are modified. 
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Recently, there has been a great deal of renewed interest in the physics of extra dimensions, the existence of which 
is a generic feature of string theories. In this paper, we consider the cosmology that results if our observable universe 
is a 3-brane in a higher dimensional universe. In particular, we focus on the case where our 3-brane is located at the 
Zi symmetry fixed plane of a Zi symmetric five-dimensional spacetime; our work applies to any 3-brane with this 
i-Q Zi symmetry. Hence this work applies to a class of models proposed by Hof ava and Witten H which compactify M 

theory on an orbifold S 1 /Z2, with one of the Eg on one orbifold fixed plane and the other Eg on the other orbifold fixed 
plane. The size of one of the extra dimensions is determined by the orbifold radius (we will call this extra dimension 
'the bulk'). In the Hof ava and Witten models, once the rest of the extra dimensions are compactified on a Calabi-Yau 
manifold with a radius smaller than the orbifold radius, the standard model (SM) fields charged under a subgroup 
of Eg, are confined to one of the two 3-branes while gravity is free to propagate in the bulk. There has been some 
effort to explore the cosmologies resulting from such a class of models which we will refer to as the Hofava- Witten 
models. Some of this work has focused on the inflationary periods or on the modifications to the particle physics 
related phenomena as a result of Kaluza-Klein excitations. There has been additional study of cosmology in scenarios 
without a Z2 symmetry. References of previous cosmological work in the context of our observable universe as a 
3-brane in a higher dimensional context include [|2|-|4|. 

In this work we focus on rather general cosmological considerations of a 3-brane with Zi symmetry. We explore 
consequences of the main novel features: i) the boundary conditions on the orbifold fixed planes induced by the Zi 
symmetry and ii) the dependence of the 4-dimensional Planck constant on the orbifold proper length. In the first 
two-thirds of this paper, we restrict ourselves to models in which the spacetime has the structure Ma x S 1 /Zi where 
M.4 is foliated by flat, homogeneous and isotropic 3-slices and find that, generically, substantial modifications to the 
observable Friedmann- Robertson- Walker (FRW) universe result. In the last third of the paper, we make remarks about 
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causality and the horizon problem that apply to any theory in which the volume of the extra dimension determines 
the four-dimensional gravitational coupling. While this work was in progress, some of the same results have been 
presented by ||; however, much of our work covers additional aspects of cosmology in theories with extra dimensions 
as wellQ 

Because the novel elements of Hof ava-Witten models result in a different set of cosmological equations, observational 
constraints can be imposed on such models. We find that, from a five dimensional point of view in general, FRW 
cosmologies are not the generic solutions on the brane containing our world. In particular, there is no manifest 
dynamical mechanism which drives the the theory to have the fields living on the brane dominating the stress energy. 
Hence, one rarely recovers the standard relationship H 2 ~ pl^v? p i where p is the energy density of the fields living 
on the brane and H is the Hubble speed on the 3-brane. The main point is that the behavior of H and thus 
the evolution of our observable universe is controlled by the modulin from the extra dimensions and the attending 
boundary conditions of the spacetime, rather than by what is on our brane. This statement is true even when the 
bulk is "empty." One can think of this as a version of the moduli problem, which states that our 3-brane's stress 
energy can be generically subdominant compared to the effect of the moduli (particularly those not confined to our 
3-brane) on our 3-brane's cosmological history. It is possible although not generic to reproduce ordinary FRW on 
our brane, and we construct a specific example in which FRW does result. Although we discuss this problem in the 
context of a compact S /Z<2 extra dimension motivated by the Hof ava-Witten model, this conclusion should apply to 
any theory in which our observable universe is a 3-brane located at the Z2 symmetry fixed point of a Zi symmetric 
five dimensional spacetime. 

The altered relationship between H and p leads to drastic modifications of the results from Big Bang Nucleosynthesis. 
The temperature/time relationship at the time of Big Bang Nucleosynthesis is drastically modified, and the element 
abundances produced in a generic model are in violation of observation. We also discuss the growth of density 
perturbations required for large scale structure formation and find that the standard results are no longer obtained: 
galaxy formation cannot proceed in the usual fashion. 

We also find the result that, in a broad class of models, it is generically impossible to stabilize the radius of the 
orbifold without causing the equation of state w = P/p on our brane to diverge; i.e., fine tuning of the equation of 
state is required in order to allow one to fix the four dimensional Planck mass today at a constant value. Because the 
3-brane observable universe that we live in has an energy density whose scale is much smaller than the fundamental 
Planck scale, we would naively expect that the stress energy of the 3-brane is irrelevant as far as the stabilization of 
the orbifold is concerned. However, this need not be true in general since the brane stress energy must be consistently 
coupled to the bulk stress energy. We give an example where an absolutely static orbifold is physically forbidden 
due to 1) this bulk/boundary coupling and 2) the fact that the equation of state on the boundary is constrained to 
satisfy — 1 < w = P/p < 1. In this class of models, we find that without fine tuning of the equation of state to 
w = —2/3 or without having just a cosmological constant type of equation of state w = —1, it is impossible to "pin" 
the four dimensional Planck mass to today's value; i.e., it is impossible to stabilize the radius of the orbifold. Instead, 
the four dimensional Planck mass must still be changing, albeit slowly, in order to satisfy the constraints and the 
equations of motion. We find that even when the Newton's constant is allowed to vary slowly (within experimental 
bounds), it is difficult to construct a realistic model that does not contain an unnaturally large number. The large 
number arises from the mismatch between the energy density on our 3-brane and the energy density set by the scale 
of the fundamental Planck's constant, i.e., a version of the cosmological constant problem. Although we have only 
shown this difficulty with stabilizing the radius of the extra dimension (and hence the four-dimensional Planck mass) 
explicitly for specific models, we suspect that this difficulty may persist more generally for other models with extra 
dimensions. 

We also make some remarks about causality and the solution to the horizon problem in theories with extra dimen- 
sions. These remarks apply to any theory, regardless of boundary conditions, in which the observable universe lies on 
a three-brane and gravity is described by an extra dimension. We caution that many assumptions go into the standard 
inflationary requirements, and illustrate how a few of these assumptions may be modified here. In particular, the 
assumption that H 2 ~ T 2 /m p i(t c ) at an early time may be modified. In addition, the fact that the value of the four 
dimensional Planck mass changes as the radius of compactification changes must not be forgotten in inflation models, 
and in fact may be exploited to solve the horizon problem (see also ||). 

In Section II, we present the model with a discussion of boundary conditions. In Section III, we discuss the resulting 
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nonstandard cosmology on our brane. In Section IV, we explore constraints from nucleosynthesis, structure formation, 
and the age of the universe in the context of the nonstandard cosmology. In Section V we show that, in a broad class 
of models, the four dimensional Planck mass cannot be fixed at a constant value. In Section VI we discuss causality 
and the horizon problem. We conclude in Section VII. 

II. MODEL 

In the first part of this paper we restrict ourselves to a scenario in which standard model particles are confined 
to one of the two orbifold fixed planes (3-branes) and gravity resides in the bulk. We consider a space-time with 
the structure S 1 /Z2 X M4 where S 1 /Z2 is an orbifold and M4 is a smooth manifold. We use the metric convention 
+,—,—,—,—. Let the action be decomposed as 

£5 = TT- 2 / d Xy/gTZ + Sorb + Sboundary + SgH (1) 

where 1Z is the Ricci scalar in 5 dimensions, g is the absolute value of the determinant of the five dimensional metric 
gMN, K5 is the five-dimensional Newton's constant («§ = l/m^ 5 ), S or b represents the orbifold (or the bulk) action, 
Sboundary represents the boundary (orbifold fixed plane) action (i.e. integrated only over the boundary), and Sgh is 
the Gibbons-Hawking boundary term to be discussed below. 

A. Metric and Einstein Equations 

We consider the metric of the form 

ds 2 = e Mr,«) dT 2 _ e 2«(r,«) dx 2 _ e W(T,u) du 2 (2) 

which foliates the space into flat, homogeneous and isotropic spatial 3-planes. Here x — x\, X2, £3 are the coordinates 
on the spatial 3-planes while u is the orbifold coordinate (this metric was considered for example by Lukas, Ovrut, 
and Waldram J|] within a similar context). Let S 1 be parameterized by u G [— 1, 1} with the endpoints identified. 
The Z 2 symmetry takes u — > —u and leaves two three dimensional planes fixed. These orbifold fixed planes then 
are at u = and at u = I. In this paper we will refer to these as the 3-branes. As in the Hofava-Witten model 
B, one of these 3-branes is our observable universe (which we choose without any loss of generality to be at u = 0) 
and the other is a hidden sector that communicates only gravitationally with our visible sector. Thus we define 
a W T )) = e Q< - r '" =0 - ) as the scale factor describing the expansion of the 3-brane that is our observable universe where 
t(r) = J drexp(^(r, u = 0)) is the proper time of a comoving observer (with the brane's position fixed at u=0; but 
see below). The four dimensional Newton's constant is inversely proportional to the distance between these 3-branes, 

Gn ~ 77 5 — ^7^r (3) 

where 

R(t) = f exp(P(T,u))du. (4) 

Jo 

Here e^ r '") is the scale factor describing the expansion of the orbifold. 

Note that we can always reparamctcrize the theory such that the 3-brane positions are fixed in orbifold coordinate 
space. However, in that case, we can no longer have the full freedom of gauge in the (t, u) plane. For example, we 
can no longer make a coordinate transformation such that v{t 1 u) — f3(r,u) (conformally flat in the r-u subspace), 
because such a transformation generally will give the 3-branes a time dependence for their position along the orbifold. 
However, since we can always reparameterize u by itself or r by itself, we can rewrite any metric having the property 
f3 — h> = i(u) or /?— v = v(t) (wherej(u) and «(r) are real functions) in a conformally flat form (in the r — u subspace) 
even with 3-brane coordinate positions fixed. 

By varying the action Eq. (fy) we find that the bulk equations are 



e- 2 "Too = 3e- 2/3 (a'/3' - a" - 2a' 2 ) + ie- 2v {a 2 + d/3) (5) 

e- 2a T u = e~ 20 {Za n - 2d '0 + 2a! V - /?V + v' 2 + 2a" + v") + 

e- 2v (-da 2 - 2d/3 ~ (i 2 + 2av + (3v - 2a - ft) 
e- 2p T M = -3(a' 2 + aVJe" 2 ^ + 3e- 2u (2a 2 - av + a) 

T 40 = 3v'a + 3a'(-a + $)-3a' (6) 

where the dots denote partial derivatives with respect to r while the primes denote partial derivatives with respect 
to u. 

An alternative and equivalent way to think about the five-dimensional space is a boundary picture. Here the 
coordinate u is restricted to one half of the circle, u € [0, 1] and five dimensional spacetime is written as [0, 1] x M.4. 
The orbifold fixed planes then turn into boundaries of this five-dimensional space. 

We will take the 3-branes to be at fixed orbifold coordinate u. However, the general results of our paper, including 
the result that ordinary FRW cosmology is modified by the existence of the bulk, would qualitatively hold even if one 
were to allow the branes to "bend", i.e., give the location of the 3-branes some w-dependence. 

B. Boundary Conditions 

The bulk equations only apply in an open region that does not include the boundary. By continuity, the metric on 
the boundary (brane) is determined by the metric in the bulk. It is the Israel conditions which connect the boundary 
and the bulk. Hence, the bulk stress energy cannot be set to zero without imposing strong constraints on the energy 
density of the boundary. 

General remarks on Israel Conditions: Here we briefly remind the reader of the Israel conditions (see [p| and 
references therein), which are relevant when there are boundaries (or equivalently discontinuities of the derivatives of 
the metric). Alternatively, one can derive equivalent conditions by just inserting a 6 function for the matter action 
(see below) as is done in [||§]. 

The Gibbons-Hawking term in Eq. (|l|) above takes the form 

Sgh = \ J2 f d 4 xVhK (7) 

K5 faces Jbux 

where K = h MN Kmn is the trace of the extrinsic curvature Kmn = h^VpriN, Hmn = gMN + timTIn is the induced 
metric on the boundary surface, and the vector n , orthonormal to the surface, points toward the region which the 
surface bounds.^] 

Upon varying S5 with respect to the metric we find, in addition to the usual 

TZ-mn — l/2<7Miv7^ = T MN (8) 

(where Tmn — -7^ s mn ) within the bulk (which was used in the last subsection), the boundary variation 

f oSboundary x MN f ^GH cMN , ' f j5„, /-jm? MN n , a \ 

J s g MN °9 + J J^mn 5 9 + 2^ J d x V9on M N9 - (9) 

where the last integrand is a total divergence and becomes a surface term. This leads to the Israel conditions 

y^ (K M n - KIimn) = t M N , (10) 

ifaces 

where the sum over faces is for each side of the boundary surface and we have defined 



3 The index M, for example, takes integer values from to 4. For example, if we work in the coordinate frame of eqn. (2) 
with the boundary at a fixed coordinate value, then the fact that n m nNg MN = 1 implies that ri4 = a/— 344 = e~ . Thus 
ft.44 = <?44 + n\ = 0, while Hmn can be nonzero for M, N < 4, so that Iimn indeed projects onto the boundary. It is then easy 
to see that only the tangential derivatives (P = 1,2,3) survive in obtaining Kmn- 
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as the energy momentum tensor on the boundary. We will assume that this energy momentum tensor on the boundary 
can be written in a perfect fluid form. For example, 

£° = 4p (12) 

and 

t\ = -4P, (13) 

where p and P are the energy density and pressure, respectively, measured by a comoving observer. One could write 
the 00 component of eqn. (10) in the following way, if one also uses eqn. (12). Defining Hi n — -& where dv = e@du 
(by analogy with the Hubble constant), we can write the 00 Israel condition as p = -^ X)±faces -^»«- 

The Israel condition in eqn. (10) is very similar to AE perp — Aira in electrostatics, where a is the charge/area on a 
charged plate and E perp is the component of the electric field perpendicular to the plate. Eqn. (10) intimately relates 
the energy-momentum on our brane to its extrinsic curvature, i.e., to the way that it is imbedded in the bulk, just as 
the charge on a plate is intimately connected to the electric field on either side. 

Note that there is an alternative way of looking at the problem which gives the same results. Instead of partitioning 
the orbifold into regions with boundaries (an approach adopted from [|3|), one could elevate the integrand of the 
surface action to a density in five dimensions by inserting a 5 function. In that case, we would obtain 

T MN = TMN + tMNS(S) (14) 

where E represents the boundary surface with the appropriate measure for the S function built in. In lieu of eqns. (8) 

and (10) we would then use IZmn — \9mnT^ = Tmn everywhere. 

Note also that the Israel conditions Eq. din) probe the global structure of the spacetime even though each of 
the conditions is local. They probe the global structure because they result from the boundary terms arising from 
integration over all space. Hence, loosely speaking, the topological constraints of the S 1 jZi are encoded into local 
boundary conditions by the Israel conditions. 

Note that no additional information, beyond what is contained in Einstein's equations and the Israel conditions, is 
obtained from considering the Bianchi identities. In the bulk, the Bianchi identities automatically give VmT mn = 0.Q 
On the boundary, V m~£ MN = X)±faces h NM TMpn p allows the energy to flow between the boundary and the bulk. 
Since these equations follow from the Einstein equations and from Eq. (JIOJ) , there are no additional constraints coming 
from these equations. 

Israel conditions for our model: We can now find the Israel conditions specific to the metric in Eq. (0). Our 
boundary conditions are found by considering the two 3-branes at u — and u — I. These are the planes at which we 
have discontinuities in the normal derivatives of the metric, and hence where the Israel conditions for the boundaries 
give nontrivial constraints. 

As mentioned above, we can always reparameterize the theory such that the 3-branes are fixed in coordinate space. 
Because we are considering a brane at a Zi symmetry fixed plane, from the sum over faces in eqn. (10) we get two 
identical terms, merely resulting in a factor of two. In a more general spacetime, the Israel conditions in eqn. (10) 
would still hold, but one would have to explicitly add the contributions from the two sides of the boundary surfaces 
since these contributions would no longer be equal. With the Z2 symmetry, the Israel conditions are 

T6a'|± = e^ 2 %o (15) 

±v'\± = \e p - 2a k l + \e^t m (16) 

where the | implies that the equation is to be evaluated at one of the two sides of the boundary surface. In each 
equation, the upper signs are to be used for the face at the boundary plane at u = + while the lower signs correspond 



4 If the bulk contains only one field (say a scalar field), then its equation of motion will be determined by the Bianchi identities; 
if there are more fields, then other equations of motion are needed as well. 



to the face at the boundary plane at u = l~ ^ Note that because of the Z-i symmetry defining the orbifold fixed planes, 
the lower sign also applies to the the surface u = Q~, for example. We have also denoted the boundary parallel spatial 
coordinates x with the indices i. Using eqn. (12) and goo = e 2l/ , we see that eqn. (15) can be written p = 6a'rnp /5 e _/3 . 
As mentioned above, if one makes a particular gauge choice such as v(t, u) = (3{t,u), then one must allow the 
possibility of a moving 3-brane. Only if f3— v is independent of time or independent of u can one choose v(t, u) = /3(r, it) 
with fixed brane positions without any loss of generality. The Israel conditions for a 3-brane having a time dependent 
coordinate r(r) are 



too = T6e 



tn = ±e 2a 



2(i/-/J) 



(e 2l V + e^ar) 



( e 2„ _ e 2/V2)VeF2? - e- 2v i 
r + 2e 2 (' y - /3 V + r 2 (/3' - 2a') + e 2{ - v '^>v' - 2fV + 2f(a + /3) - e 2 ^-^r 3 (/3 + 2d) - rv 

( e 2v _ e 2{Sj,2\^/ e -2f3 _ e -2vj.2 



(17) 
(18) 



As explained before, with a proper choice of coordinates, the 3-brane positions can be fixed (r = 0) at the coordinates 
that we have chosen for all bulk coordinate time. However, if we insist instead that v = (3, then this apparent loss of 
a functional degree of freedom will reappear in the nonzero r except under certain conditions previously discussed. 

III. NONSTANDARD EVOLUTION OF OUR 3-BRANE 

^From the 4-dimensional point of view of a low energy observer living on the 3-brane, the effect of the extra 
dimensions can be considered as the coupling of some extra moduli fields to the boundary energy density. The 
distinguishing feature of the brane scenario is that these extra fields arise from the gravitational coupling to the 



extra dimensions and hence must obey the boundary equations Eq. (15) and Eq. (|16|). From Eq. (]5|), not only is the 
usual four-dimensional Friedmann equation relating the energy density and the Hubble speed modified, but Eq. ( jig ) 
restricts the possible relationship even further. In general, there seem to be no restrictions of the stress energy tensor 
that drive the energy density confined on the 3-brane to behave as H 2 (an explicit illustrative example of this will be 
given later in the next section V). Hence, the relationship between the energy density of matter confined to the brane 
and the expansion scale factor will generally be different from the standard cosmology. 

Let us consider why one does not generically recover the standard evolution for our observable universe: i.e. why in 
general H 2 will not be proportional to the energy density confined to the observable sector 3-brane. The main point 
is that the behavior of H and thus the time evolution of our observable universe can be controlled by the moduli from 
the extra dimensions and the attending boundary conditions of the spacetime, rather than by what is on our brane. 
One can think of this as a version of the moduli problem, which states that our 3-brane's stress energy is subdominant 
compared to the effect of the moduli coming from extra dimensions on our universe's history. 

The above statement is most transparent when one considers the full five dimensional Einstein's equations rather 
than the dimensionally reduced 4-dimcnsional effective action. For example, consider the case where there is no 
stress-energy in the bulk, i.e. T^ v = 0. Still, the behavior of the extra dimension is determined by the 3-brane 
stress energy t^ Vl and the derivatives along the extra dimensions of metric components act as energy sources for 
the Friedmann equation governing the expansion of the observable 3-brane. Consider, for example, Eq. (|5|). From 
a higher dimensional point of view, the first of the equations in Eq. (O) clearly shows the contributions to H that 
are independent of Too, the bulk energy density. For example, we can take /3 to be a constant ($ = (3' = 0) and 
set Too = a" = 0. Then, using the Israel condition in Eq. (|l5|), as well as Eq. (||), Eq. (|l2|), and the fact that 
H = a < = e~ v a, we obtain the result 

H= 8 -^L(V2lp), (19) 

which shows a linear relationship between p and H . Hence, even if the bulk is empty, we generally will not recover 
the standard FRW cosmology. This unusual relationship has also been discussed by Rcf . H , which considered the 4-4 
component of the Einstein's equation with the assumption 



3 Here, for example, u = + corresponds to approaching u — from the positive side (u > 0). 



Pbuik « %=*, (20) 

where pbrane is the energy density on the brane and pbuik is the energy density on the bulk, and showed that an 
unusual form of the Friedmann equation ensues (a linear relationship between H and p). 

In general, almost any relationship between p and H is possible (e.g. such as p ~ H q ). One can see this, for example, 
by using the fact that a" and j3' on the 3-brane are not well constrained. From the four dimensional effective theory 
point of view, this is the same as saying that the moduli not confined to our brane can contribute to the energy 
density such that nearly arbitrary time dependence for the expansion scale factor a may be achieved. More explicitly, 
suppose we seek a solution of the form 

P = pH q (21) 

where p is a constant which we can parameterize as M 4: ~ q c where M is a mass scale and c is a dimensionless constant. 
In that case, the Israel equations, at u = for example, become 

- i/(r, u = 0) = 3a' (r, u = 0)(w(t) + |) (22) 

6a' (t, u = 0)e~ /3(T ' u=0) = ^«§(e- y(T ' u=0) d(r, u = 0)) q (23) 

where u> = — i^/^o is the equation of state. Note that because these constraint equations are only evaluated at a point 
on the orbifold, one can consider v{t,u = 0) and u'(t,u — 0) as independent functions. Similarly, a(r,u = 0) and 
a'(r, u — 0) can be considered independent functions. 

For example, we can construct a solution for any value of q with w = P j p = (pressureless fluid) and p ex 
l/a 3 (r, u = 0) = exp(— 3a(r, u = 0)) as follows: 



0{t,v,)=v(t,u)='-£ ^-Y pu (24) 

L3rJ 
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3r 
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a(r )U ) = — ^ f- (F(u)-F(« = 0))+ln(A) + ^M-) (25) 



3 r ' 



where -F(u) is any smooth function satisfying F'(u = 0) = 1 and A, r , and p are constants. For example, suppose 
q = 2 (corresponding to the usual FRW relationship) and F(u) — u. The resulting stress energy tensor components 
are 



-4exp(^gT)(-8lT 4 + 4^ 2 (t 2 + u 2 )) 
243r 6 



T o _ —J^PK 27rJ )K-° rl -r_^5P V t» I) , 0R , 

1 o — KT^Zs \ > 



Tl -4exp(^r) K I M (-36r 2 n + K lMr 2 -4 U 2 )) 

Jl _ 24376 [27) 

-2exp(^gT)(81r 4 + 54 K 2 M r 2 M + 2^ M 2 (t 2 + 16 M 2 )) 
T4 2437^ (28) 

T4 ° - ^437^ (29) 

Ap 

9r 2 

P = (31) 



P=S ( 3 °) 



Hence, with a nontrivial time dependent bulk stress energy tensor in the bulk, one can recover the FRW expansion 
on the boundary. Note that in this particular example, there is no energy flowing between the brane and the bulk. 

Of course, once a specific model is written down with specific choices for the stress-energy tensor components T^^, 
we will be able to solve Einstein's equations together with the boundary conditions to calculate exactly which q results 
in p ~ H q . At present we can only point out that different choices of T M „ will give rise to different values of q, and 
generically will produce modifications to the standard FRW results. 

Let us now consider the modifications to the standard FRW cosmology from a purely four dimensional effective 
action point of view. The total action, Eq. (fil) , can be dimensionally reduced to 

£4 ~ 7p 2 / ^ x V~94{'R"l + AT?.) + S or bi + Sboundaryl + Shoundaryl (32) 



where A1Z is the gravitational moduli contribution that arises from dimensionally reducing the 5-dimensional Ricci 
scalar, S or b4 is the contribution coming from the bulk field, and the remaining action terms correspond to the 
brane-confined field contributions. In the Hofava-Witten model, Sboundaryi corresponds to the contribution from the 
observable sector and Sboundary2 corresponds to the contribution from the hidden sector. What we are calling the 
moduli problem is that Sboundary2, Sorbi, and ATZ contributions can dominate the energy density determining the 
expansion rate of the universe. For a particular orbifold geometry about which the dimensional reduction is taken, the 
stress energy contributions coming from ATZ, Sboundaryi, Sboundary2 and Sorbi are related through the Israel conditions 
and the Einstein equations. In Ref. |5J], a linear relationship between H and p arising from this moduli dominance 
was emphasized, particularly when Sorbi — 0. We would like to merely point out that in general, any relationship 
between H and p may result if all the moduli contributions (i.e. ATZ, Sboundaryii and Sorbi) are taken into account. 

IV. CONSTRAINTS 

In the last section, we saw that the ordinary matter confined to our brane does not (alone) determine the expansion 
rate of our universe H; in contrast the dilution of the energy density and the pressure of the ordinary matter proceeds 
as usual as the universe expands. Hence, any physical observable that measures the expansion rate independently of 
the "ordinary" energy density will be able to test this 3-brane scenario. To illustrate this constraint, we will consider 
big bang nucleosynthesis (BBN), structure formation, and the age of the universe. 

Big Bang Nucleosynthesis: As discussed above, generically one obtains modifications to the cosmological standard 
model relationship between evolution of the scale factor and energy density on our brane. These modifications will 
drastically alter nucleosynthesis. In BBN, a given reaction rate depends directly on the energy density (the photon 
number density, for example) while the freeze out is governed by the ratio of the Hubble speed to the reaction rate. 
Hence, the energy density time variation can be compared with the universe expansion rate since the ratios of the 
various light element abundances measure the various temperatures at which different elements freeze out. 

At the high temperatures in the early universe, the ratio of neutrons to protons is determined by its thermal 
equilibrium value, 

n/p = e- Q/T ,T>T F , (33) 

where the neutron-proton mass difference Q = 1.293 MeV. Neutrons drop out of equilibrium below a freeze-out 
temperature Tp, where the weak interaction rates can no longer keep up with the expansion of the universe. Below 
Tf the n/p ratio continues to fall due to /3-decay on the time scale of the neutron half-life r„. In the standard BBN 
scenario, nucleosynthesis begins at a temperature approximately given by 

MeV , , 

T D = 2.2—- 34 

— 111!) 

where r\ is the baryon to photon ratio. Once Tjj is reached, deuterium becomes stable against photodissociation and 
nucleosynthesis takes place very rapidly, efficiently converting essentially all of the available neutrons into He 4 . In 
this approximation, the primordial helium abundance Y p is given by 

2n 2n 2e~ rtD 

Yr> = ( ) n = ( ) „exp\-T(t D - t F )] f— r (35) 

p y n+p' D K n+p' F P[ y n l + cxp[Q/T F ] V ; 

where the final approximation is valid since T _1 = T„/ln2 3> t F ~ lsec. 

Now we consider the modification of these results when the equation H 2 ~ pjmz^ is modified. As an example, let 
us consider the following modification to standard 4-dimensional Einstein equations: 

H = p/M 3 , (36) 

where M is some mass scale, perhaps generically the Planck mass. Such a relationship was discussed above and is 
also found in the toy model discussed in the next section. It was also found by Lukas, Ovrut, and Waldram Q in 
their nonlinear case as well as by Q in a similar context. If we assume that the scale factor expansion is a power law, 
we then have H ~ 1/t. In a radiation dominated universe, p ~ T 4 . With Eq. (pq), one can then find that the scale 
factor grows as a ~ i 1 / 4 and the modified temperature-time relation is T ~ i~ 1 ' ¥ M 3 / 4 . Recall that the freeze-out of 
neutrons and protons occurs when a typical n *-* p weak interaction rate T ~ GpTp is equal to the expansion rate 
H. Here H ~ T 4 /M 3 (as opposed to the usual T 2 /m p i of the standard BBN model). Thus we find 



T F = T F (m pl /M 3 ) 1 (37) 

where an overbar indicates the standard BBN model value. For example, if M = m p i — 10 19 GeV, we have Tp = 
10~ 44 Tf- This model has the novel feature that Td > Tf, so that the exponential factor in the numerator of Eq. ( |35| ) 
is not there. However, both temperatures are significantly lower than the usual Tp, so that the production of He^s 
driven to zero! 

In addition to the specific case considered above of H oc p, one can allow almost any q in the equation p oc H q , 
depending on the T^ v in the bulk. To analyze a variety of possibilities, we consider power law growth of the scale 
factor both in the matter and the radiation dominated regimes, and allow the power law to be a free parameter. Thus 
we take the scale factor on the 3-brane to be e a<yT,u= ^ ~ (t(r, u — 0)/ti) n for the radiation dominated era, where U 
is a constant and t(r,u — 0) is the proper time of the comoving observer on the 3-brane. During matter domination 
we take the power law index to be to instead of n (with to and n usually different). Consider the neutron-proton 
interconversion interactions freezing out. By setting the reaction rate equal to the Hubble expansion rate, to within 
an order of magnitude, the freeze out temperature then is found to be 



(rpl/m \ (l/m-5) 

Gp^(po m /por) [1/n - 1/m] j (38) 

where To is the temperature of the cosmic background photons today H$ is the Hubble speed today, po m is the matter 
density today, and po r is the radiation density today. 

Compare for example the two possible cases: 1. {n = 1/2, m — 2/3} (usual FRW scenario) and 2. {n = 1/4, m = 
1/3} (a possible brane scenario). For case 1, we obtain a crude estimate Tp ~ 1 MeV while for case 2, we obtain a 
crude estimate of Tp ~ 10 14 GeV. Hence, even if this crude estimate were orders of magnitude off,[] we will obtain a 
large measurable discrepancy for the resulting element abundances since the neutron to proton ratio is constrained 
to a precision of less than 1%. 

Structure Formation: Consider the Jeans analysis of a fluid overdensity within a universe with a ~ t m during the 
matter dominated epoch. To obtain an idea of the type of effects that we may find, instead of analyzing the exact 
perturbation equations with metric perturbations includedjj we will merely modify the time dependence of the scale 
factor in the usual Jeans analysis equation valid for standard FRW cosmology. We assume the usual long wavelength 
limit and consider only one component nonrelativistic matter fluid. Normalizing the background energy density and 
the Hubble speed to approximately what we measure today (i.e. p « 3/(8tt)MZHq), the fluid overdensity equation 
can be written as 

S"(x) + —5'{x) - ^x- 3m 6 = (39) 

x 2 

where 5 = (p — p)/p is the fluid overdensity and x = t/to with to denoting the time today. This equation can generally 
be solved in terms of Bessel functions for constant to. In the standard FRW cosmology with matter domination, 
to = 2/3, and there is one solution to S which grows as the scale factor S ~ a ~ t 2 / 3 and the other decreasing as 1/x. 
However, take for example the case where p ~ H which is equivalent to setting to = 1/3 if p oc a~ 3 . In this case, 
there will be two growing solutions for <5, one having the form 1 + x/A + ... and the other a; 1/ ' 3 (l + x/8 + ...). This 
agrees with the expectation that in a universe that is expanding more slowly, the overdensity will grow more quickly 
with the scale factor. 

Extra growth of the perturbations compared to the standard picture might enable some new models to be considered, 
that are ruled out in the standard model. For example, one might be able to better tolerate a smooth component 
in the universe (which ordinarily suppresses perturbation growth) and still get enough growth of the fluctuations 
to get structure today. For example, particles that decay to radiation might be able to lead to a longer radiation 
dominated universe subsequent to the initial matter/radiation equality. However, the standard picture of cold dark 
matter might well be in trouble with altered growth of the scale factor. With this altered growth, calculation of the 
integrated Sachs Wolfe effect and normalization to the COBE data would undoubtedly lead to anisotropies in the 
power spectrum which are less scale invariant than the standard model, such that one would not be able to match 
the data on all different scales. 



6 The actual freeze out temperature will be lower, making the proton-neutron relative abundance ratio (which is approximately 
exp(— lMeV/Tp)) discrepancy from the observed smaller. 
7 This more careful treatment is in progress. 



The time dependence of the scale factor for a given type of ordinary energy density p can be tested by considering 
the location of where the collapse and diffusion processed power spectrum starts to flatten out. This flattening occurs 
in the standard CDM scenarios because a particular nonrelativistic matter fluctuation Fourier mode grows only after 
the matter domination phase begins and the mode wavelength enters the horizon. The horizon length at the time of 
matter radiation equality is approximately t eq . For wavelengths longer than t eq (at that time), the modes are outside 
of the horizon and will not grow until they enter the horizon. Hence, today, the mode's amplitude will be suppressed 
with respect to the mode that grew starting at time t eq . Modes that had already entered the horizon by the time t eq 
will grow together, and hence the spectrum will be flat in that region. 

With a different time dependence for the scale factor a ~ t m , the location of this bend in the power spectrum will 
be shifted from Ai to A2 

A a = ^f^)-" I A 1 . (40) 

I \P0MJ 

For the case p ex H (m=l/3), we have A2 = 10~ 7 Ai if we assume that pom is order of the critical density (i.e., 
8.1 x l(T 47 GeV 4 ). 

Age of the Universe With the time dependence of the scale factor altered, the age of the universe may be modified. 
For example, if we assume that structures form during a period after the matter-radiation equality with a ~ t m and 
that the scale factor at matter-radiation equality aeq is much smaller than today's scale factor ao, i.e. (ao/aeq) 3> 1, 
then the time period of structure formation is given by t = m/H where H is the Hubble speed today. Hence the 
lower bound on the age of the oldest stars in a globular clusters may be used to constrain m. Taking the value of |9[| 
of t > 10.2 Gyr, we find 

m > h (41) 

where the Hubble speed is H — lOOh Km/s/Mpc. This would rule out p ex H, since in that case m = 1/3, whereas 
observations of the Hubble constant find that h > 0.5. 

V. DIFFICULTIES OF STATIONARY PLANCK MASS COSMOLOGY 

We have discussed the altered relationship between the energy density and the Hubble speed in the 3-brane scenario 
compared to the standard cosmology. Another difference from the standard cosmology is that Newton's gravitational 
coupling constant can have a time dependence in the Hofava-Witten scenario because it is inversely proportional to 
the orbifold proper length. Naively, one would, in constructing a cosmological scenario, then try to freeze the orbifold 
radius. However, as we will see, this turns out to be unphysical, at least in the case where the metric can be written 
in a conformally flat form in the t — u subspacc. If one allows the orbifold length to vary slowly but within the limits 
of astrophysical bounds, we still expect some constraint although it is not as severe. We will explicitly illustrate this 
with a toy model. This difficulty of freezing the four dimensional Planck's constant is yet another manifestation of 
the cosmological constant problem. 

We will consider the stress-energy tensor in the comoving frame. As in Eq. (n2h and Eq. (|13|), we denote the energy 
density and the pressure as t[j = K i/° an d t\ = —k§P respectively. Hence, Eq. ( Jlq ) can be rewritten as 

v'\ ± = ±(\w+ l -)ty (42) 

where we have introduced a time dependent function w = P j p defining the equation of state for the stress energy. 
Now, suppose we have a metric with the property that v — /3 is purely a function of u or r; any such metric can be put 
into the conformally flat form in the r — u subspace. In other words, for any such metric we can effectively take v = /3 
in Eq.(2). Given such a metric, we will now attempt to fix the orbifold radius and consequently the four-dimensional 
Planck mass by taking fi to be time independent. Since Eq. (M3) holds for all time, we can impose the condition 

j/ l±= ° (43) 

where, since /3 is time independent, t = e^"V. Solving Eq. (f43| ) and Eq. ( ff2[ ) for w{t), we find 

-2 w(U) + 2/3 

w{t) = ~ + IpWpW) ' (44) 
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where ti represents some initial time. 

Hence, for all energy densities that decrease with time and w(U) ^ —2/3, our equation of state grows in magnitude 
without bound. For the usual macroscopic stress energy consisting of homogeneous scalar fields and fermionic particles 
bound to the brane, the equation of state is constrained to — 1 < w(t) < 1. Hence, Eq. (J44|) most likely will violate 
any realistic equation of state, and the four dimensional Planck mass cannot be fixed to a constant value (e.g. with 
a potential). We comment briefly on the possibility that w(t) — w(ti) — —2/3 exactly, the only value that allows 
the radius to stabilize for decreasing energy density. We note that this value has been proposed (Wang et al [|10[) 
as the best fit value to the combination of a number of data sets, including recent supernovae observations as well 
as various measures of large scale structure and the microwave background. However, in the context of Eq. (0), 
obtaining this value would require either some unknown dynamical mechanism to drive the number to w = —2/3, 
or an uncomfortable fine tuning. We can see that, for a universe with a decreasing energy density, the denominator 
in the second term continues to decrease with time and drives the second term away from the phenomenologically 
acceptable value —1 < w < 1, unless w(ti) in the numerator is fine tuned to be arbitrarily close to —2/3. 

We wish to comment that the difficulty we find in attempting to stabilize the radius of the extra dimensions to a 
fixed value is similar to the problem of trying to have a static Einstein universe. Both require fine tuning. While a 
static universe can be achieved with a proper parameter choice, it is unstable to perturbations about these values; 
similarly, a stabilized radius of the extra dimension is unstable to slight deviation of w(ti) from the precise value of 
2/3. Although we have only been able to prove the fine tuning required for radius stabilization for specific metrics in 
the Hofava Witten scenario, we suspect that such a result may generalize to any attempt to stabilize the radius in 
the context of an extra dimension in which gravitons propagaten. 

A way to avoid these conclusions and still to satisfy Eq. ( p4| ) is if p(t) on the observable 3-brane is dominated by 
a constant vacuum energy (an effective cosmological constant) for all time. Then p(t) = p(ti), w(t) — w(ti) = — 1, 
and Eq. (M4) is trivially satisfied. Such an energy density on the brane does not automatically lead to exponential 
expansion of the brane if there is a nonlocal contribution canceling the effect of the 3-brane cosmological constant, 
and hence this case can possibly be made to match today's universe. However, generically Eq. ffl) disagrees with the 
evolution of our observable universe. 

The condition that the metric on the orbifold be absolutely static is clearly too stringent. The radius of the orbifold 
must be allowed to change. Then the only experimental constraint on the variation of the orbifold length comes 
from the bound on the variation of the Newton's constant (see for example |L2 ] and references therein) . The most 
conservative bound is approximately 

dG N ,n__ ^ a-., m-ll™-l 



dt 



/G N < 6 x lO^yr , (45) 



which comes from the spin down rate of pulsars. Although this bound is only measured on astrophysical time scales, 
we will assume it to be valid for cosmological time scales. We remind the reader that Gn is inversely proportional 
to the distance between 3-branes (Eq. (0)), which is obtained via Eq. M) as an integral over the entire bulk. Hence 
the best way to avoid the constraints of Eq. (J43|) and Eq. (J44J), while still satisfying the observations on dGjq /dt, is to 
have a large time variation of the orbifold metric near the brane and vanishing time variation elsewhere throughout 
the bulk. We will construct such an example shortly. 

The most general constraint from Eq. (|16| ) arises as a result of w being of order unity. This implies that, at u = 0, 

\e-P0l\ ~ e = \t° l\ ~ l(T m (46) 

where I ~ 10~ 15 GeV -1 is the physical length scale of the orbifold appropriate for the Hofava- Witten model considered 



in 14]. The numerical estimate in Eq. (46) is made by taking the energy density to be p ~ 10~ 29 gm cm~ a as can be 
roughly inferred from rotation curves of galaxies (we motivate this number with rotation curves rather than with the 
typical critical density of standard cosmology because the latter may change in these models). This condition says 
that, very near the brane, the orbifold scale factor remains approximately independent of the transverse distance off 
the brane. Below we will allow the physical length scale of the orbifold / to be as large as a mm and find that our 
conclusions are relatively unchanged. 



Now we will construct an explicit example which demonstrates that the time variation bound and Eq. (46) can be 
simultaneously satisfied by allowing the orbifold length to have a time variation. However, as we will see, the model 
is still unnatural. Consider the following ansatz: 



Remarks on the difficulty of stabilizing the radius in extra dimensions can also be found in pi] , which appeared shortly after 
this work 

if 



P=^e- M/2 cosh[\{u-l/2)] (47) 

A 

where \l ^> 1. Note that this ansatz fits the optimal situation described below Eq. fl45| ) in that most of the time 
variation is near the boundary: dGiy/dt is determined by dfl/dt with /? (and its time derivative) biggest at the 
boundaries. We define the orbifold proper length to be 

l 
e 0(r(t),u) du (48) 



where t is the comoving time on the brane. Then, from equations Eq. (0), Eq. (Eq), Eq. (|47J), and Eq. (Eq), the time 
variation bound on the Newton's constant translates into 



where C\ — 6 x fO n yr 1 . Dropping factors of order unity, combining Eq. (49) and Eq. (|46|), we find the constraint 
on h to be 

2 < \h\ (50) 



c^nn 



where e = lO" 114 for Z - l(r 15 GeV _1 . For example, if one takes t = 10 10 yr, then if \h\ > 10 233 GeV/yr, the time 
variation constraint on the Newton's constant is satisfied. Note that if one took I as large as 1 mm, the parameter 
will still be unnatural; since eocl, the bound on h ex l/l 3 and we would require \h\ > 10 149 GeV/yr. Of course, such 
a large dimensionful constant most likely is unphysical. Incidentally, for I ~ 10~ 15 GeV _ , A > 10 127 GeV which is 
much greater than l/l as was assumed below Eq. (|47]). Note also that although Eq. ( (49J ) seems to imply an upper 
bound on the value of \h\, instead a lower bound arises because A in Eq. (jig) has to adjust as \h\ is varied to keep 
the equation of state reasonable. This requirement of such an unphysically large number for h is a result of t^l being 
very small. This in turn is just a restatement of the cosmological constant problem in the sense that there is a severe 
mismatch between the natural scale 1/(k 2 1) and the energy density on the brane. One way of accomodating the severe 
mismatch is to have the long age of the universe (corresponding to the small energy density of the universe) built into 
the solution itself. Indeed, Eqs. (E4) and ( |25| ) with q = 1 give such an example, where 1/r ~ O(l/10 10 yr) supplies 
the unnaturally small (compared to 1/(k 2 Z) t ' 4 ) energy scale. 

A cosmological example of our 3-brane for v = {3 

Now we turn to evolution of the 3-brane when v — (3. We consider the other Israel equation Eq. (jl5|) which specifies 
p once a is given. Suppose we have the brane expanding as a power law in comoving time. Thus, we write 

a = lnf(u) + mln(t(T,u)/U) (51) 

where m is the exponent of the power law and ti is a constant. Using the ansatz for j3 in Eq. ([47J), we then find 

p = A(u)jln(l+te) + B(u) Y ^ § (52) 

where A(u) and B{u) are some functions of u and 

6 = 0/t. (53) 

Note that, as t — > oo, p — ► (A + B/ff)/t oc H . This is in contrast with the FRW cosmology where p oc 1/i 2 ex H 2 
independently of the equation of state, as long as the scale factor is a power law with respect to time. Hence, we 
have a novel feature of the energy density scaling like the Hubble speed instead of its squared. As discussed in the 
previous section, such a relationship between the energy density and the time causes many cosmological problems. 
Furthermore, if one writes down the equation of state corresponding to this cosmology, one will see that w approaches 
a constant value that is determined, e.g., by the energy density at some initial time. 

Note that this behavior relies strongly on the assumption v = /3, which cannot always be made to be true if the 
orbifold fixed planes have time independent coordinate positions. Indeed, without the choice v = (3, the physical time 
derivative taken to obtain the Hubble speed has no general relationship with the physical spatial (orbifold directional) 
derivative taken to obtain the energy density p. Note that the bulk equations are largely unconstrained because of 
the freedom we have in choosing a(r, u) for u away from the boundaries. Even with the simple power law expansion 
of the scale factor discussed in this section, f(u) in Eq. ( plj ) is unconstrained except at u = 0. 
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VI. CAUSALITY AND THE HORIZON PROBLEM 

Interesting new avenues for solutions to the horizon problem arise in the context of extra dimensions. Here we 
confine ourselves to a few remarks and reserve a more complete discussion to a work in progress. Previous work by 
Freese and Levin Jl3| addressed the horizon problem in some generality, of course in a four-dimensional context. We 
discuss here one aspect of the generalization of that work to higher dimensions. The major point we wish to make is 
that the application of standard inflation requirements, such as the requirement of 60 e-foldings, may be modified in 
the context of extra dimensions as some of the assumptions that went into obtaining this requirement are relaxed. 

Note that the discussion in this section is true for any model in which the length of the extra dimension(s) determines 
the four-dimensional gravitational coupling. For example, it applies to periodic boundary conditions considered by 
many authors including fl4j| . 

An estimate of the size of the observable universe today is given by the distance light could travel between photon 
decoupling and now, d i, s ~ a(t ) L ° dt/a(t) = 0(1) X (t — tdec) for a oa t p and p = 0(1) between tdec & n d t . Here 
a is the scale factor. We can compare the comoving size of the observable universe to the comoving size of a causally 
connected region at some earlier time t c : dhor(t c )/a(t c ) — L c -4|y. The observable universe today fits inside a causally 
connected region at t c if 

tthorKpc) ^ &obs (ZA\ 

a(t c ) ~ a(t ) 

Here, subscript o refers to today and subscript c to some early time. If condition Eq. rt54) is met, then the horizon 
size at t c (before nucleosynthesis) is large enough to allow for a causal explanation of the smoothness of the universe 
today. Note that more creative explanations of large scale smoothness may not involve comparing these two patches. 
For example, in the context of the brane scenarios, one might imagine that two regions of our observable universe 
which seem to be causally disconnected might in fact have talked to each other because of a geodesic between them 
that went off our brane, into the bulk, and then back onto our brane at some distant point. Such alternatives are 
certainly worth investigating further. In the remainder of our discussion here we restrict ourselves to the case where 
Eq. ( p3 ) is relevant; this is certainly the case for all the brane and boundary inflation models considered to date. 
For power law expansion of the scale factor both before t c and after tdec (which may or may not be the case), we 



can take t ~ H during these periods. The causality condition Eq. (54) then becomes 



1 >^r- (55) 



a c H L 

We take the Hubble constant today to be 

H = al^/mpiito), (56) 

where a = (8tt /3)(tt 2 /30)g*(t o )r] o , g* is the number of relativistic degrees of freedom and rj(t ) ~ 10 4 — 10 5 is the 
ratio today of the energy density in matter to that in radiation. 

Following Freese and Levin p3], one can further simplify this condition by making a series of simplifying assumptions. 
First, we assume for now that the entropy on the 3-brane both at time t c and today obeys S ex (aT) 3 such that we 
can write a in Eq. (|54|) in terms of the entropy. A work in progress is reexamining this condition. 

A further assumption is always made in obtaining requirements for inflation to succeed. It is that the Hubble speed 
at the early time t c is given by 

H c ~T?/m pl {t c ), (57) 

where T c is some mass scale, possibly the temperature, at time t c . As can be seen by the earlier results in this 
current paper, this assumption may not always hold. Thus one should in principle always go back to the original 
equation Eq. ( p4[ ) to check that any particular inflation proposal really does solve the horizon problem. However, for 
the models that have been heretofore proposed, the usual requirements used are at least qualitatively similar to the 
actual requirements implied by above equations. 

Even if above assumptions hold, yet another fact must be considered: the changing four dimensional Planck mass. 
Using Eqs. (pq), (pq), (p7[), and the entropy relation, the causality condition becomes 

^MA)i/B> 10 - 2 ^. (58) 

m p i(t ) A c Jo 
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One can now invoke the idea of changing the four dimensional Planck mass in lieu of changing the entropy in order to 
solve the horizon problem, as has been suggested by Levin and Freese |15| and as is studied in the context of string 
theory by Veneziano [[16| and collaborators. Here, in the context of extra dimensions, if the compactified dimension 
that determines the strength of gravity shrinks with time, one may be able to use the changing four dimensional 
Planck mass to explain causality as in Eq. ( |5q ) (see also Riotto J8|). 

In the context of inflation, one must also be careful about the changing four dimensional Planck mass due to the 
changing size of the extra dimension in obtaining the appropriate causality condition. Suppose inflation ends at time 
t e . Assuming that entropy increases by the end of inflation while it remains constant after inflation, we can then 
rewrite the condition Eq. ( pq ) as 

^>10- 2 -^(^4). (59) 

a c l e l x m p i(t c )' 

It is this last factor of m p i(t ) / rn p i(t c ) that must not be forgotten, if the value of the four dimensional Planck mass 
changes between the time of inflation and now. One way to think about this is that in taking our observable universe 
today and rescaling it to some early time (to see if it fits inside a causal radius), the size of the observable patch may 
be bigger (or smaller) than one would naively think because of the difference in the Planck masses. For example, in 
the inflation model of |L4| inflation takes place at an early stage while the radius of compactification is very small and 
the Planck mass much smaller than today. If one takes m p i(t c ) =TeV, for example, then one needs the scale factor 
to grow by an additional factor of 10 19 GeV/TeV = 10 16 compared to the standard inflationary result with constant 
Planck mass. Instead of the usual 60 e-foldings one needs 100. This is not a significant difference, since it is usually 
just as easy to get 100 e-foldings as 60. Still, we caution that there are many assumptions that go into the calculation 
of requirements for inflationary models, and these should in principle be checked. 

VII. CONCLUSIONS 

In this work, we have explored the general cosmology of a flat, homogeneous and isotropic 3-brane located at the 
Zi symmetry fixed plane of a Zi symmetric five dimensional spacetime as in the Hofava-Witten model compactified 
on a Calabi-Yau manifold. We have found an important manifestation of the moduli problem: the usual relationship 
between the FRW scale factor and the energy density on our 3-brane is modified from p ~ H 2 to almost any relationship 
between p and H, such as p ~ H q with q ^ 2. Hence, in these models, one does not generically recover the standard 
FRW evolution for our observable universe. The time evolution of our observable universe can be controlled by the 
moduli from the extra dimensions and the attending boundary conditions of the spacetime, rather than by the energy 
density on our brane. These modifications persist even in the limit that the bulk stress energy vanishes, because of the 
modifications to Einstein's equations in the bulk and the nontrivial Israel conditions arising from the Zi symmetry. We 
have shown that this typically leads to difficulties with primordial nucleosynthesis, structure formation, and bounds 
on the age of the universe. 

We have also found that, in a large class of cases, it is impossible to stabilize the radius of the extra dimension 
(i.e., "pin" the 4-dimensional Planck mass at a fixed value) without fine tuning the equation of state on our brane. 
The difficulty we find in attempting to stabilize the radius of the extra dimensions to a fixed value is similar to the 
problem of trying to have a static Einstein universe. We find that the radius of the extra dimension can be frozen 
at a fixed value in a universe (on the 3-brane) with decreasing energy density only if the equation of state satisfies 
w = P I p = —2/3; however, the radius is unstable to slight deviations from this value such that w must be fine-tuned to 
extreme precision to exactly this number. Although this particular value has recently been discussed, e.g. by Wang et 
al jlQ] as tentatively favored by overlap of a number of astrophysical observations, the lack of a dynamical mechanism 
to drive w to this value renders this exact number extremely fine tuned. Even a slowly varying Planck mass, in 
agreement with observations, requires a tuning. This tuning is a result of the energy density of the universe today 
being severely mismatched with the energy density set by the scale of the fundamental Planck constant. Although we 
have only been able to prove the fine tuning required for radius stabilization for specific metrics in the Hofava-Witten 
scenario (those with conformally flat metric in the t — u subspace), we suspect that such a result may generalize to 
any attempt to stabilize the radius in the context of an extra dimension in which gravitons propagate. 

We repeat, however, that this constraint involving the equation of state can be evaded if one assumes that the 
energy density confined to the 3-brane, p, is dominated by the cosmological constant. In that case, the equation of 
state as well as the energy density can remain constant in time. A cosmological constant which dominates the energy 
density p confined to the brane does not necessarily lead to an inflationary universe; in principle there can be other 
(nonlocal) contributions from off the brane canceling the effects of the cosmological constant, such that the evolution 
of our universe today can proceed as usual. The ordinary matter in this case would contribute only negligibly to p. 
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We can apply the results of our paper to the model of Randall and Sundrum |17|], which falls into the class of 
models we have considered here. First, our result that FRW cosmology does not generically result on our brane would 
apply in their model also. Second, they presumably avoid the problems of stabilizing the orbifold radius by having a 
cosmological constant dominate the energy density of our brane, as described in the previous paragraph. 

Finally, we made some remarks about causality and the horizon problem in theories with extra dimensions. We 
discussed possible modifications to the standard inflationary picture. For example, the time variation of the gravita- 
tional coupling due the changing size of the compactified dimension can be used to solve the horizon problem or at 
least modify the number of e-foldings needed to solve the horizon problem within the context of inflation. In addition, 
the modifications indicated here to the standard FRW relation H 2 oc p would change the requirements of inflationary 
models. 
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